Electrically tunable g-factors in quantum dot molecular spin states 
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We present a magneto-photoluminescence study of individual vertically stacked InAs/GaAs quan- 
tum dot pairs separated by thin tunnel barriers. As an applied electric field tunes the relative ener- 
gies of the two dots, we observe a strong resonant increase or decrease in the g-factors of different 
spin states that have molecular wavefunctions distributed over both quantum dots. We propose a 
phenomenological model for the change in g-factor based on resonant changes in the amplitude of 
the wavefunction in the barrier due to the formation of bonding and antibonding orbitals. 
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Quantum Dots and Quantum Dot Molecules (QDMs) 
have proven to be a versatile medium for isolating and 
manipulating spins P, Q, which are of great interest 
for quantum information processing 0, 0|. In partic- 
ular, photoluminescence (PL) spectra have been used 
in self-assembled QDMs to observe coherent tunneling 
0, 0, li and identify spin interactions through fine 
structure |9( . Electrical control of isolated spins through 
their g-factors is highly desireable for implementation of 
quantum gate operations. To date, electrical control of 
g-factors has only been observed in ensembles of electrons 
in quantum wells by shifting the electron wavefunctions 
into different materials [lO|, IllL Il2l Il3| . In this Letter we 
present a striking electric field resonance in the g-factor 
for molecular spin states confined to a single quantum 
dot molecule. 

To our knowledge this is the first observation of elec- 
trical control over the g-factor for a single confined spin. 
Moreover, the isolation of a single QDM allows us to 
spectrally resolve and identify individual molecular spin 
states that have different g-factor behaviors. In Fig. ^ 
we indicate molecular spin states of both the neutral ex- 
citon (X", one electron recombining with one hole) and 
positive trion {X~^, electron-hole recombination in the 
presence of an extra hole) at zero magnetic field. The 
different electric field dependences of the g-factors for 
these states is apparent in Fig.^, where the splitting of 
PL lines increases for some molecular spin states and de- 
creases for others. This electric field dependence is nearly 
an order of ma gnit ude lar ger than previously reported in 
quantum wells [lO|, IllL \l2L . The effect arises from the 
formation of bonding and antibonding orbitals, which re- 
sults in a change in the amplitude of the wavefunction in 
the barrier at resonance. 

Our QDMs consist of two vertically stacked self- 
assembled InAs dots truncated at a thickness of 2.5 nm 
and separated by a 2 nm GaAs tunneling barrier [l^ . 
As an applied electric field tunes the relative energies of 
the two dots, strong tunnel coupling between the hole 
states creates the molecular spin states. Unlike samples 
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FIG. 1: (a) B = T photoluminescence spectra of a single 
QDM. The complex pattern of anticrossings arises from the 
formation of molecular spin states, (b) At B = 6 T, the molec- 
ular spin states have a Zeeman splitting (bars) that depends 
strongly on the applied electric field (F). 



with a thicker tunnel barrier , the states retain molec- 
ular character throughout the observed range of electric 
fields. We present data from a single molecule, but the 
universality of the behavior has been verified by detailed 
studies of 7 other molecules from the same sample. We 
first explain the spectra and molecular spin states at B 
= T. We then describe the magnetic field dependence 
and propose a phenomenological model for the electric 
field dependent Zeeman splitting. 

In Fig. 121 we show all PL lines from X° and X+ at B = 
T. These lines are identified by their relative energies, 
the power dependence of their intensities and the electric 
field dependence of the anticrossings j^- The X° lines 
(Fig. |3i,) show a clear anticrossing at Fxo. The anti- 
crossing arises from tunnel coupling between the ground 
state hole levels in each dot, which forms bonding (Xg) 
and antibonding (AT^) orbitals. The electron remains in 
the bottom dot throughout the range of electric fields 
considered here H. The X+ lines (Fig.|3)) have a more 
complicated pattern because anticrossings occur in both 
the initial (one electron and two holes) and final (one 
hole) states 

To explain the origin of the A+ molecular spin states 
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FIG. 2: Energies extracted from Fig. (a) lines anti- 
cross at Fxo, where the direct (lower inset) and indirect (up- 
per inset) transitions are degenerate, (b) initial states 
anticross at Fx+ and near 90 kV/cm. Final (hole) states an- 
ticross at F^. Inset: charge distribution for the circled singlet 
and triplet transitions. 



we turn to the Hamiltonians. The basis states will be 
identified as /j^'/j^-^^, where es [ex] are the spins of 
electrons (±1/2: f or |) and Hb [hr] the spins of heavy 
holes (±3/2: -ff or JJ.) in the bottom [top] dot. X indicates 
an exciton (h a single hole) and Q is the net charge, k is 
the total spin projection. Singlets, which have total spin 
±1/2, will be denoted Xs to distinguish them from the 
±1/2 triplets. 

The final state has only a single hole with relative en- 
ergies given by the diagonalization of the Hamiltonian: 




(1) 



in the spin degenerate ^'o''-+3/2' o'i|'''+3/2 basis, where t}^ 
is the tunneling energy for a bare hole, d is the distance 
between dot centers and F is the applied electric field. 
The energies of the final states as a function of electric 
field are shown in Fig.lSJ). The formation of bonding and 
antibonding states at the anticrossing point is illustrated. 

The initial state (X+) contains one electron and two 
holes. For simplicity we present only the electron-spin- 
up case, which is degenerate with the spin-down case at 
zero magnetic field. If both holes are in the same dot, the 
Pauli principle requires singlet configurations: |||,'o^5 



(both in the bottom dot) or Q^^j^Xg (both in the top dot). 
If there is one hole in each dot, singlet A-'^^g ) and triplet 



(l'l^+i/2' i,4^-5/2' l;°^^7/2) configurations are pos- 
sible. By |''J|.^5 we mean the antisymmetric hole spin 



X^, ,„ means the sym- 



wavefunction (-d-jft" — tf,!^). 
metric hole spin wavefunction (J|, lt + lt, J|). The relative 
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FIG. 3: Calculated zero magnetic field energies of spin- 
degenerate (a) initial states (X"*") and (b) final states {h'^) 
[13 ■ Labels indicate the dominant basis state for easy com- 
parison to Eq. 121 Insets in (b) show the bonding (lower) and 
antibonding (upper) wavefunctions (green) and probability 
distributions (red), (c) Fine structure of singlet and triplet 
states at F = 79.2 kV/cni. Arrows indicate optically allowed 
transitions to the bonding final state. 



energies of the initial states are given by the Hamiltonian: 
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M.O^^S' ^.t^^-S' 4,t"^ + l/2' ^.4"^-5/2' ^,-^^ + 7/2' 

t-ll-^'^'S basis. Fi and F2 are due to Coulomb interac- 
tions, tx+ is the Coulomb-corrected tunneling energy 
for a hole in the presence of an electron and addi- 
tional hole, Ek = edF + mJ'^'^ with m = 0,-1-1,-1 for 
k = -fl/2,-5/2,+7/2. J"^ is the exchange e nerg y be- 
tween an electron and hole in the same dot Di- 
agonalizing Eq.[2| gives the energies of the initial states, 
which are plotted as functions of electric field in Fig. Oi. 

Because tunneling is a spin conserving process, only 
the Itxt singlet state can tunnel couple with |^'g^5 



and 



T, 



xi 



P^jj^^vg, which must be singlets because the two 
holes are in the same dot. These 3 singlets are therefore 
strongly mixed to create molecular orbital states that an- 
ticross near 57 kV/cm {Fx+) and 90 kV/cm. Unlike the 
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FIG. 4: (a) Magnetic field dependence of initial (X"*") and 
final (h) states for the singlet and triplet transitions at an 
electric field of 79.2 kV/cm (schematic Fig.j^t). Vertical lines 
indicate spin allowed recombinations, (b) Calculated (colored 
lines) and experimentally observed (black points) PL energies, 
(c) Zeeman splitting of PL lines (A) and initial and final states 
for ^^5/2 (red lines) and Xg (green lines). 



singlet states, the triplet states do not mix and are not af- 
fected by these anticrossings ^3 ■ This creates a "kinetic" 
splitting between triplet and singlet states Til . An exam- 
ple is indicated by the dashed oval in Fig. [33, where the 
lower energy singlet line remains separated from the two 
optically allowed triplet lines, which are split by electron- 
hole exchange. The fine structure of the corresponding 
states is shown in Fig. 

By adding the Zeeman interaction to the Hamiltonian 
[l9j |. we calculate the magnetic field dependence of the 
molecular spin states. Fig. ^ shows the states of Fig. 

which split into doublets with an applied longitudinal 
magnetic field. The final states are simply the two spin 
orientations of a single hole, split by g^. The splitting of 
initial states depends on their spin configuration. Due to 
the two parallel hole spins, X^^^^ has a large splitting 

given by ge + Igt while ^±7/2 '^^ ^P^^^ ~5e + "^gh- 
In contrast, the ±1/2 singlet {Xg) and triplet [X'^-^^^) 
have oppositely paired hole spins and therefore a small 
splitting given by g,,- The g-factor for PL transitions is 
given by the difference in g-factor between the initial and 
final states. Away from the electric field resonances, the 
g-factor is gT = ge + dh, as indicated by the vertical lines. 

Using the model described below, we obtain gx = 
—2.2. To plot the initial and final states in Fig. 0Ji we 
have taken the relative weights of the electron and hole 
g-factors to match those obtained by Bayer: ge ~ —0.6 
and gii = — 1.6 0|. The calculated energies of the tran- 
sitions are shown by the lines in Fig. ^p. The experi- 
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FIG. 5: Zeeman splitting and corresponding g-factor, gr, as 
a function of electric field at B = 6 T. 

mentally observed PL energies are given by the symbols, 
with the diamagnetic shift (10.9 /xeV/T^) subtracted. 

The g-factor resonances are clearly evident in Fig. [SJ 
where the symbols plot the measured energy splitting of 
the X'^ and two X+ Zeeman doublets at B = 6 T as a 
function of electric field. Strong enhancement or suppres- 
sion of the splitting is observed at Fji and F^o . All of the 
data can be qualitatively explained by a phenomenolog- 
ical model of the formation of bonding and antibonding 
orbitals, which results in resonant changes in the ampli- 
tude of the wavefunction in the barrier. We focus first on 
transitions involving a bonding orbital (open symbols), 
which have an increased amplitude of the wavefunction 
in the barrier (Fig. [SJa lower inset). 

The wavefunction for the bonding orbital of a single 
hole can be written as "^b = a|l) + b\2), where |1) and 
2) are the basis states of Eq. ^ the wavefunctions for 
holes localized in the two different dots. The coefficients 
a and b are functions of electric field determined by Eq.^ 
The electric field-dependent g-factor for a hole in a bond- 
ing orbital is given by g^{F) — B\gh{z)\^ b) , where 
gh{z) is the hole g-factor as a function of position in the 
sample. gh{z) is taken as a phenomenological parame- 
ter, in part because the degree of alloying between the 
nominally pure InAs dots and GaAs barrier is unknown 

If we assume that the g-factors for the holes localized 
in each dot are the same, we get g^{F) = g^ + 2abgi2, 
where gi2 = (l|f//i(z)|2) gives the contribution from the 
amplitude of the wavefunction in the barrier. In the case 
of the Xg singlet (shown in Figs. and 2];), the initial 
states are split by g^, so the total transition g-factor is 
given by 

gB^F)=ge+g^{F)=gT+ . (3) 

Je^d^iF-Fn)^+Ul 
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where the explicit form for 2ab determines the hneshape 
centered around the anticrossing point Fh- 

The lower black line in Fig.Elis obtained by fitting Eq. 
13 to ^Xg. Using the measured values of (0.86 meV) 
and Fji (63.3 kV/cm) we find gi2 = 1.65. The barrier 
contribution is positive, like the heavy hole g-factor in 
bulk GaAs (~ 1.05) |2^|23|. Because gr and gi2 have 
opposite sign, the splitting reaches a minimum at F^, 
where the amplitude in the barrier is at a maximum. 
The g-factor at the minimum is —0.44. 

Because the ±5/2 triplet states recombine to the same 
bonding orbital of the final hole (Fig. 0), the model pre- 
dicts the splitting of these lines as a function of elec- 
tric field with no additional fitting. As shown in Fig. 
0];, the ±5/2 triplet states have an initial state splitting 
of ge + 2gh. The transition g-factor is therefore given 
by griP) = (ffe + '^9h) - {gh + "^abg^) = gr - 2a6gi2. 
This is shown by the upper black line, which matches 
AX^g^2: the observed splitting of the ±5/2 triplets. The 
maximum splitting corresponds to a g-factor of —4.23. 

To apply the g-factor model to the bonding orbital of 
the neutral exciton (lower branch Fig. [2Jt) we use the 
measured tunneling energy (tjfo=0.58 meV) and anti- 
crossing field (Fxo=82.1 kV/cm). The lower red line in 
Fig. El shows the fit to AX'^, the Zeeman splitting of 
PL lines from the bonding orbital. We find g[2 = 1-32. 
The electron-hole Coulomb interaction is responsible for 
the difference in tunneling energy and anticrossing field 
from the bare hole case and is also likely responsible for 
the difference between gi2 and g[2- The g-factor at the 
minimum is —0.59. 

Using this value of g'i2-, the model immediately explains 
the increase in splitting for the antibonding orbital (up- 
per branch Fig. Et-); which has a reduced wavefunction 
amplitude within the barrier. The g-factor for the anti- 
bonding orbital is given by g^{F) = gx — 2a5(7^2: which 
increases in magnitude at the resonant field because gr 
is negative and 175.2 positive. This is shown by the up- 
per red line in Fig.|3 which matches AX^, the Zeeman 
splitting for excitonic recombination from the antibond- 
ing orbital. The splitting increases to a maximum (g- 
factor —3.35) at the anticrossing point. The antibonding 
transitions for X+ show similar behavior, but arc too 
weak to obtain full resonance curves. 

The overall agreement between the model and exper- 
imental data is quite good. There are some minor dis- 
crepancies, which highlight the need for a detailed the- 
ory, possibly requiring inclusion of excited states |23l |. 
However, the agreement of the data with the resonance 
linewidths calculated using independently measured val- 
ues of txo and th is strong confirmation that the g-factor 
dependence does arise from the formation of bonding 
and antibonding orbitals. For X+ the g-factor resonance 
arises from the wavefunction of the single hole in the final 
state, while for the orbital wavefunction of the hole 
is influenced by the additional electron. 



We also studied samples in which electrons tunnel 
through the barrier, with an anticrossing energy (~ 2.3 
meV) comparable to that of the hole tunneling sample 
presented here (~ 1.7 meV). This requires a thicker bar- 
rier (10 nm) because of the smaller electron effective 
mass. The electron wavefunction amplitudes in the bar- 
rier should be at least as large as the hole-tunneling case. 
However, the electron g-factor in bulk GaAs (—0.44) 
is similar to the electron g-factor in InAs quantum dots 
(—0.6) and according to our model the contribution from 
the barrier should not significantly change the electron 
g-factor. We see no electric field dependence of the g- 
factor in these electron anticrossing samples. By adding 
aluminum to the barrier, the resonant contribution to the 
electron g-factor could be enhanced. 

We have presented a resonant change in g-factor as a 
function of electric field for the molecular spin states of 
QDMs with a thin tunnel barrier. By studying single 
QDMs, we are able to identify the individual molecular 
spin states and the different resonant behavior of their 
g-factors. The results suggest that design of molecular 
spin states and tunnel resonances may provide new op- 
portunities for combining optical and electrical control of 
confined spins. 
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